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We study the Fermi level structure of 2 + 1-dimensional strongly interacting electron systems in 
external magnetic field using the AdS/CFT correspondence. The gravity dual of a finite density 
fermion system is a Dirac field in the background of the dyonic AdS-Reissner-Nordstrom black 
hole. In the probe limit the magnetic system can be reduced to the non-magnetic one, with 
Landau-quantized momenta and rescaled thermodynamical variables. We find that at strong 
enough magnetic fields, the Fermi surface vanishes and the quasiparticle is lost either through 
a crossover to conformal regime or through a phase transition to an unstable Fermi surface. In 
the latter case, the vanishing Fermi velocity at the critical magnetic field triggers the non-Fermi 
liquid regime with unstable quasiparticles and a change in transport properties of the system. We 
associate it with a metal-" strange metal" phase transition. Next we compute compute the DC Hall 
and longitudinal conductivities using the gravity-dressed fermion propagators. For dual fermions 
with a large charge, many different Fermi surfaces contribute and the Hall conductivity is quantized 
as expected for integer Quantum Hall Effect (QHE). At strong magnetic fields, as additional Fermi 
surfaces open up, new plateaus typical for the fractional QHE appear. The somewhat irregular 
pattern in the length of fractional QHE plateaus resemble the outcomes of experiments on thin 
graphite in a strong magnetic field. Finally, motivated by the absence of the sign problem in 
holography, we suggest a lattice approach to the AdS calculations of finite density systems. 
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I. INTRODUCTION 

The study of strongly interacting fermionic systems at finite density and temperature is a challenging task in 
condensed matter and high energy physics. Analytical methods are limited or not available for strongly coupled 
systems, and numerical simulation of fermions at finite density breaks down because of the sign problem There 
has been an increased activity in describing finite density fermionic matter by a gravity dual using the holographic 
AdS/CFT correspondence [I]. The gravitational solution dual to the finite chemical potential system is the electrically 
charged AdS-Reissner-Nordstrom (RN) black hole, which provides a background where only the metric and Maxwell 
fields are nontrivial and all matter fields vanish. In the classical gravity limit, the decoupling of the Einstein-Maxwell 
sector holds and leads to universal results, which is an appealing feature of applied holography. Indeed, the celebrated 
result for the ratio of the shear viscosity over the entropy density Q is identical for many strongly interacting theories 
and has been considered a robust prediction of the AdS/CFT correspondence. 

However, an extremal black hole alone is not enough to describe finite density systems as it does not source the 
matter fields. In holography, at leading order, the Fermi surfaces are not evident in the gravitational geometry, but 
can only be detected by external probes; either probe D-branes [I] or probe bulk fermions [3-01- Here we shall consider 
the latter option, where the free Dirac field in the bulk carries a finite charge density Q. We ignore electromagnetic 
and gravitational backreaction of the charged fermions on the bulk spacetime geometry (probe approximation). At 
large temperatures, T ^ fi, this approach provides a reliable hydrodynamic description of transport at a quantum 
criticality (in the vicinity of superfiuid-insulator transition) At small temperatures, T fi, in some cases sharp 
Fermi surfaces emerge with either conventional Fermi-liquid scaling [sf or of a non-Fermi liquid type Q with scaling 
properties that differ significantly from those predicted by the Landau Fermi liquid theory. The non-trivial scaling 
behavior of these non-Fermi liquids has been studied semi-analytically in and is of great interest as high- Tc 
superconductors and metals near the critical point are believed to represent non-Fermi liquids. 
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What we shall study is the effects of magnetic field on the holographic fermions. A magnetic field is a probe 
of finite density matter at low temperatures, where the Landau level physics reveals the Fermi level structure. The 
gravity dual system is described by a AdS dyonic black hole with electric and magnetic charges Q and H, respectively, 



esp e 

corresponding to a 2 + 1-dimensional field theory at finite chemical potential in an external magnetic field [10[ . Probe 
fermions in the background of the dyonic black hole have been considered in [ill [l^ ; and probe bosons in the same 
background have been studied in [l3|. Quantum magnetism is considered in jl4| . 

The Landau quantization of momenta due to the magnetic field found there, shows again that the AdS/CFT 
correspondence has a powerful capacity to unveil that certain quantum properties known from quantum gases have 
a much more ubiquitous status than could be anticipated theoretically. A first highlight is the demonstration (l5j 
that the Fermi surface of the Fermi gas extends way beyond the realms of its perturbative extension in the form of 
the Fermi-liquid. In AdS/CFT it appears to be gravitationally encoded in the matching along the scaling direction 
between the 'bare' Dirac waves falling in from the 'UV boundary, and the true IR excitations living near the black 
hole horizon. This IR physics can insist on the disappearance of the quasiparticle but, if so, this 'critical Fermi-liquid' 
is still organized 'around' a Fermi surface. The Landau quantization, the organization of quantum gaseous matter 
in quantized energy bands (Landau levels) in a system of two space dimensions pierced by a magnetic field oriented 



in the orthogonal spatial direction, is a second such quantum gas property. We shall describe here following [11 
that despite the strong interactions in the system, the holographic computation reveals the same strict Landau-level 
quantization. Arguably, it is the mean- field nature imposed by large N limit inherent in AdS/CFT that explains this. 
The system is effectively non-interacting to first order in 1/A^. The Landau quantization is not manifest from the 
geometry, but as we show this statement is straightforwardly encoded in the symmetry correspondences associated 
with the conformal compactification of AdS on its flat boundary (i. e., in the UV CFT). 

An interesting novel feature in strongly coupled systems arises from the fact that the background geometry is only 
sensitive to the total energy density + contained in the electric and magnetic fields sourccd by the dyonic black 
hole. Dialing up the magnetic field is effectively similar to a process where the dyonic black hole loses its electric 
charge. At the same time, the fermionic probe with charge q is essentially only sensitive to the Coulomb interaction 
gqQ. As shown in flTj . one can therefore map a magnetic to a non- magnetic system with rescaled parameters (chemical 
potential, fermion charge) and same symmetries and equations of motion, as long as the Reissner-Nordstrom geometry 
is kept. 

Translated to more experiment-compatible language, the above magnetic-electric mapping means that the spectral 
functions at nonzero magnetic field h are identical to the spectral function aX h = Q for a reduced value of the coupling 
constant (fermion charge) q, provided the probe fermion is in a Landau level eigenstate. A striking consequence is that 
the spectrum shows conformal invariance for arbitrarily high magnetic fields, as long as the system is at negligible to 
zero density. Specifically, a detailed analysis of the fermion spectral functions reveals that at strong magnetic fields 
the Fermi level structure changes qualitatively. There exists a critical magnetic field at which the Fermi velocity 
vanishes. Ignoring the Landau level quantization, we show that this corresponds to an effective tuning of the system 
from a regular Fermi liquid phase with linear dispersion and stable quasiparticles to a non-Fermi liquid liquid with 
fractional power law dispersion and unstable excitations. This phenomenon can be interpreted as a transition from 
metallic phase to a "strange metal" at the critical magnetic field and corresponds to the change of the infrared 
conformal dimension from > 1/2 to < 1/2 while the Fermi momentum stays nonzero and the Fermi surface 
survives. Increasing the magnetic field further, this transition is followed by a "strange-metal" -conformal crossover 
and eventually, for very strong fields, the system always has near-conformal behavior where kp = {) and the Fermi 
surface disappears. 

For some Fermi surfaces, this surprising metal-" strange metal" transition is not physically relevant as the system 
prefers to directly enter the conformal phase. Whether a fine tuned system exists that does show a quantum critical 
phase transition from a FL to a non-FL is determined by a Diophantine equation for the Landau quantized Fermi 
momentum as a function of the magnetic field. Perhaps these are connected to the magnetically driven phase transition 
found in AdS5/CFT4 [l^. We leave this subject for further work. 

Overall, the findings of Landau quantization and "discharge" of the Fermi surface are in line with the expectations: 
both phenomena have been found in a vast array of systems [l3l and are almost tautologically tied to the notion of 
a Fermi surface in a magnetic field. Thus we regard them also as a sanity check of the whole bottom-up approach of 
fermionic AdS/CFT [i-i,!!!, giving further credit to the holographic Fermi surfaces as having to do with the real 
world. 

Next we use the information of magnetic effects the Fermi surfaces extracted from holography to calculate the 
quantum Hall and longitudinal conductivities. Generally speaking, it is difficult to calculate conductivity holograph- 
ically beyond the Einstein-Maxwell sector, and extract the contribution of holographic fermions. In the semiclassical 
approximation, one-loop corrections in the bulk setup involving charged fermions have been calculated [isj . In an- 
other approach, the backreaction of charged fermions on the gravity-Maxwell sector has been taken into account and 
incorporated in calculations of the electric conductivity [l]. We calculate the one- loop contribution on the CFT side. 
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which is equivalent to the holographic one-loop calculations as long as vertex corrections do not modify physical 
dependencies of interest [l^, [l^ . As we dial the magnetic field, the Hall plateau transition happens when the Fermi 
surface moves through a Landau level. One can think of a difference between the Fermi energy and the energy of the 
Landau level as a gap, which vanishes at the transition point and the 2 + 1-dimensional theory becomes scale invariant. 
In the holographic D3-D7 brane model of the quantum Hall effect, plateau transition occurs as D-branes move through 
one another In the same model, a dissipation process has been observed as D-branes fall through the horizon 
of the black hole geometry, that is associated with the quantum Hall insulator transition. In the holographic fermion 
liquid setting, dissipation is present through interaction of fcrmions with the horizon of the black hole. We have also 
used the analysis of the conductivities to learn more about the metal-strange metal phase transition as well as the 
crossover back to the conformal regime at high magnetic fields. 

We conclude with the remark that the findings summarized above are in fact somewhat puzzling when contrasted to 
the conventional picture of quantum Hall physics. It is usually stated that the quantum Hall effect requires three key 
ingredients: Landau quantization, quenched disorder ^ and (spatial) boundaries, i. e., a finite-sized sample [2^. The 
first brings about the quantization of conductivity, the second prevents the states from spilling between the Landau 
levels ensuring the existence of a gap and the last one in fact allows the charge transport to happen (as it is the 
boundary states that actually conduct). In our model, only the first condition is satisfied. The second is put by hand 
by assuming that the gap is automatically preserved, i. e. that there is no mixing between the Landau levels. There 
is, however, no physical explanation as to how the boundary states are implicitly taken into account by AdS/CFT. 

The paper is organized as follows. We outline the holographic setting of the dyonic black hole geometry and 
bulk fermions in the section |IT1 In section IIIII we prove the conservation of conformal symmetry in the presence of 
the magnetic fields. Section IIVI is devoted to the holographic fermion liquid, where we obtain the Landau level 
quantization, followed by a detailed study of the Fermi surface properties at zero temperature in section [V] We 
calculate the DC conductivities in section I VII and compare the results with available data in graphene. In section 
IVIIl we show that the fermion sign problem is absent in the holographic setting, therefore allowing lattice simulations 
of finite density matter in principle. 

II. HOLOGRAPHIC FERMIONS IN A DYONIC BLACK HOLE 

We first describe the holographic setup with the dyonic black hole, and the dynamics of Dirac fermions in this 
background. In this paper, we exclusively work in the probe limit, i. e., in the limit of large fermion charge q. 

A. Dyonic black hole 

We consider the gravity dual of 3-dimensional conformal field theory (CFT) with global U{1) symmetry. At finite 
charge density and in the presence of magnetic field, the system can be described by a dyonic black hole in 4- 
dimensional anti-de Sitter space-time, AdS4, with the current in the CFT mapped to a U{1) gauge field Am in 
AdS. We use ii,v,p, . . . for the spacetime indices in the CFT and M, N, . . . for the global spacetime indices in AdS. 

The action for a vector field A]\j coupled to AdS^ gravity can be written as 

s.^^l^.^,{n + ±-^n,„F"'-). (1) 

where gp is an effective dimensionless gauge coupling and R is the curvature radius of AdS4 . The equations of motion 
following from eq. (jlj arc solved by the geometry corresponding to a dyonic black hole, having both electric and 
magnetic charge: 

ds' = gMNdx^'dx^ = ^ {-fdt^ + dx^ + dy^) + (2) 
The rcdshift factor / and the vector field Am reflect the fact that the system is at a finite charge density and in an 



^ Quenched disorder means that the dynamics of the impurities is "frozen", i. e. they can be regarded as having infinite mass. When 
coupled to the Fermi hquid, they ensure that below some scale the system behaves as if consisting of non-interacting quasiparticles only. 
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external magnetic field: 

02 + ii-2 j^^ 



1 



„4 „3 



At = ^ (l - ^) , = hx, = Ar = 0, (3) 

where Q and H are the electric and magnetic charge of the black hole, respectively. Here we chose the Landau gauge; 
the black hole chemical potential /.i and the magnetic field h are given by 

with ro is the horizon radius determined by the largest positive root of the redshift factor /(ro) = 0: 

Af^r3 + ^!±£!. (5) 

The boundary of the AdS is reached for r — )■ oo. The geometry described by eqs. describes the boundary theory 

at finite density, i. e., a system in a charged medium at the chemical potential [i — /ibh and in transverse magnetic 
field h = /ibh, with charge, energy, and entropy densities given, respectively, by 

-0-9— ((,\ 

The temperature of the system is identified with the Hawking temperature of the black hole, Th ^ |/'(ro)|/47r, 

Since Q and H have dimensions of [i]^, it is convenient to parametrize them as 

g' = 3r^, Q2^i?2 = 3^4^. (8) 
In terms of tq, ?'* and r.^,^ the above expressions become 

J - ' W 

with 



^^= Sgp-^, h^Vigp — -■ (10) 

/t^ro R* 

The expressions for the charge, energy and entropy densities, as well as for the temperature are simplified as 

2^/3 r2 _ Ir^ + Sr^Ao _ 27r 

T-l^Al-'-k)- (11) 



47r i?2 



'0 



In the zero temperature limit, i. e., for an extremal black hole, we have 

r = ^ ro=r„, (12) 

which in the original variables reads + = ir^. In the zero temperature limit (|12p . the redshift factor / as given 
by eq. ([9]) develops a double zero at the horizon: 

/-6 ("~:"^' +0((r-n.)^) . (13) 



5 



As a result, near the horizon the AdS/i metric reduces to AdS2 x with the curvature radius of AdS2 given by 

i?2 = -^R. (14) 

This is a very important property of the metric, which considerably simplifies the calculations, in particular in the 
magnetic field. 

In order to scale away the AdS^ radius R and the horizon radius tq, we introduce dimensionless variables 

r — > r* — > ror*, r** — > ror**, 

M^r^M, Q^rlQ, ,H^rlH, (15) 

and 

[t,x)^ — (t,x), Am^-^Am, u}~^—u}, 
^^IT^^"^ ^^R-^^^ ^^W-^^ 

ds^ R^ds^. (16) 

Note that the scaling factors in the above equation that describes the quantities of the boundary field theory involve 
the curvature radius of ^^54, not AdS2- 
In the new variables we have 



At = /i ^1 - , fi = Vigprl = gpQ, h = gpH, (17) 
and the metric is given by 
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ds' = {-fdt^ + dx^ + dy^) + ^^, (18) 

with the horizon at r = 1, and the conformal boundary at r — > oo. 

At r = 0, r** becomes unity, and the redshift factor develops the double zero near the horizon. 

As mentioned before, due to this fact the metric near the horizon reduces to ^c?S'2 xR^ where the analytical calculations 
are possible for small frequencies However, in the chiral limit m = 0, analytical calculations are also possible in 
the bulk AdS4 (2l| . which we utilize in this paper. 



B. Holographic fermions 

To include the bulk fermions, we consider a spinor field "0 in the ^^54 of charge q and mass m, which is dual to an 
operator O in the boundary CFT^ of charge q and dimension 

A = ^ + niR, (20) 

with mR > — 5 and in dimensionless units corresponds to A = | + m. In the black hole geometry, eq. ([2]), the 
quadratic action for tp reads as 

S^=i J d^xy/^ (i/ir^^I?MV' - mi/ji;) , (21) 

where i/' = V'^r*! 8'iid 

■Dm =dM + ^ WahMr"^ - iqAM , (22) 
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where uJabM is the spin connection, and F"'' = i[r'',r'']. Here, M and a,b denote the bulk space-time and tangent 
space indices respectively, while fx, v are indices along the boundary directions, i. e. M = {r, /i). Gamma matrix basis 
(Minkowski signature) is given by eq. (jAl2|) as in 

We will be interested in spectra and response functions of the boundary fermions in the presence of magnetic field. 
This requires solving the Dirac equation in the bulk 0, @ : 

{t'^'Vm - m) = 0. (23) 

From the solution of the Dirac equation at small lu, an analytic expression for the retarded fermion Green's function 
of the boundary GFT at zero magnetic field has been obtained in Q ■ Near the Fermi surface it reads as : 

Gi.{^.k) = ^T"^";), (24) 

uj — vpK± — Zj(a;, I ) 

where fc_L = A: — fc^^ is the perpendicular distance from the Fermi surface in momentum space, hi and vp are real 
constants calculated below, and the self-energy S = Si + iT,2 is given by Q 

J:{uj,T)/vp = T^'^gi^) = (27rr)^-fe2e-^— - ^ ]; ^ ^ (25) 

-'^ ( 2 ^ 2irT + 6 ) 

where v is the zero temperature conformal dimension at the Fermi momentum, v = v^p, given by eq. ([58]). = f-iq, 
/i2 is a positive constant and the phase 9 is such that the poles of the Green's function are located in the lower half 
of the complex frequency plane. These poles correspond to quasinormal modes of the Dirac Equation (|23p and they 
can be found numerically solving F{uj^) ~ [l^l, with 

Fi.) . ■ k,e—i2.Tr 

^V2^'^ 27rT ^ 6 / -^V2 2ttT ^ 6 ' 

The solution gives the full motion of the quasinormal poles a;*"''(fcj^) in the complex lo plane as a function of k^. It 
has been found in 0, [22|, that, if the charge of the fermion is large enough compared to its mass, the pole closest to 
the real ui axis bounces off the axis at fcj^ = (and w = 0). Such behavior is identified with the existence of the Fermi 
momentum kp indicative of an underlying strongly coupled Fermi surface. 

At T = 0, the self-energy becomes T"^^ g{ijj /T) — > Cfco;^'^, and the Green's function obtained from the solution to the 
Dirac equation reads f\ 



Gnin, k) = , ^"^^^^ ■ , , (27) 



where k± = ^/k^ — kp. The last term is determined by the IR AdS2 physics near the horizon. Other terms are 
determined by the UV physics of the AdSi bulk. 

The solutions to p3)) have been studied in detail in [B-0I- Here we simply summarize the novel aspects due to the 
background magnetic field (formal details can be found in the Appendix). 



The background magnetic field h introduces a discretization of the momentum (see Appendix A for details): 



k fccff = ^/2\qh\l, with I e N, (28) 

with Landau level index I [l^, . These discrete values of k are the analogue of the well-known Landau levels 
that occur in magnetic systems. 

• There exists a (non-invertible) mapping on the level of Green's functions, from the magnetic system to the 
non-magnetic one by sending 



iH,Q,q)^ I^O,VQ^TH^,q^l--^^^j . (29) 

The Green's functions in a magnetic system are thus equivalent to those in the absence of magnetic fields. To 
better appreciate that, we reformulate eq. p9p in terms of the boundary quantities: 



(h,,,q,T)^ [0,^,q,T[l-^]], (30) 
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where we used dimensionless variables defined in eqs. (|15ll7p . The magnetic field thus eflFectively decreases 
the coupling constant q and increases the chemical potential fJ- = gpQ such that the combination fig = fiq is 
preserved . This is an important point as the equations of motion actually only depend on this combination 
and not on fi and q separately [ll|. In other words, eq. ([50)1 implies that the additional scale brought about by 
the magnetic field can be understood as changing fi and T independently in the effective non-magnetic system 
instead of only tuning the ratio /i/T. This point is important when considering the thermodynamics. 

• The discrete momentum k^s = y^2\qh\l must be held fixed in the transformation (j29p . The bulk-boundary 
relation is particularly simple in this case, as the Landau levels can readily be seen in the bulk solution, only to 
remain identical in the boundary theory. 

• Similar to the non- magnetic system , the IR physics is controlled by the near horizon AdS2 x geometry, 
which indicates the existence of an IR CFT, characterized by operators Oi, I €z N with operator dimensions 
,5=1/2 + 1//: 

1 / / .2\qh\l\ ^ 

in dimensionless notation, and fiq = fiq. At T = 0, when r^^ = 1, it becomes 

= ^6 {m^ + 2\qh\l)~fil. (32) 

The Green's function for these operators Oi is found to be ^/'(w) ~ w^"' and the exponents i^i determines the 
dispersion properties of the quasiparticle excitations. For ly > 1/2 the system has a stable quasiparticle and a 
linear dispersion, whereas for v < 1/2 one has a non- Fermi liquid with power-law dispersion and an unstable 
quasiparticle. 



III. MAGNETIC FIELDS AND CONFORMAL INVARIANCE 



Despite the fact that a magnetic field introduces a scale, in the absence of a chemical potential, all spectral functions 
are essentially still determined by conformal symmetry. To show this, we need to establish certain properties of the 
near-horizon geometry of a Reissner-Nordstrom black hole. This leads to the AdS2 perspective that was developed 
in 0. The result relies on the conformal algebra and its relation to the magnetic group, from the viewpoint of the 
infrared CFT that was studied in 0] ■ Later on we will see that the insensitivity to the magnetic field also carries over 
to AdSi and the UV CFT in some respects. To simplify the derivations, we consider the case T = 0. 



A. The near-horizon limit and Dirac equation in AdS2 

It was established in that an electrically charged extremal AdS'-Reissner-Nordstrom black hole has an AdS2 
throat in the inner bulk region. This conclusion carries over to the magnetic case with some minor differences. We 
will now give a quick derivation of the AdS2 formalism for a dyonic black hole, referring the reader to 0] for more 
details (that remain largely unchanged in the magnetic field). 

Near the horizon r = r,, of the black hole described by the metric ([2]), the redshift factor /(r) develops a double 
zero: 

/(r)=6^^— ^+0((r-r..)')- (33) 



Now consider the scaling limit 



c 

In this limit, the metric ([2]) and the gauge field reduce to 



732 

A-^, t = A-V, A ^ with r, C finite. (34) 



ds^ = ^{-dr^+dC') + ^^{dx'+dy') 

Ar = A^^Hx (35) 



8 



where R2 = The geometry described by this metric is indeed AdS2 x M^. Physically, the scaling limit given in 
cq. p4p with finite r corresponds to the long time limit of the original time coordinate t, which translates to the low 
frequency limit of the boundary theory: 

^ 0, (36) 

where w is the frequency conjugate to t. (One can think of A as being the frequency uj). Near the AdSi horizon, we 
expect the AdS2 region of an extremal dyonic black hole to have a CFTi dual. We refer to for an account of this 
AdS2/CFTi duality. The horizon of AdS2 region is at C — 00 (coefficient in front of dr vanishes at the horizon in cq. 
((35)) ) and the infrared CFT (IR CFT) lives at the AdS2 boundary at C = 0. The scaling picture given by eqs. p4ll35p 
suggests that in the low frequency limit, the 2-dimensional boundary theory is described by this IR CFT (which is 
a CFTi). The Green's function for the operator O in the boundary theory is obtained through a small frequency 
expansion and a matching procedure between the two different regions (inner and outer) along the radial direction, 
and can be expressed through the Green's function of the IR CFT Q. 

The explicit form for the Dirac equation (|A28p in the magnetic field is of little interest for the analytical results 
that follow; for completeness, we give it in the Appendix. Of primary interest is its limit in the IR region with metric 
given by cq. ([35]) : 



^.^dc -m+ -^a^ L + ^) - -^^^ ] ^^(0 . 0, 
/9CC V-9tt \ J y/gTitcr^M, 



(37) 



where the effective momentum of the l-th. Landau level is A; = y^2\qh\l, fj.q — iiq and we omit the index of the spinor 

field. To obtain eq. (|37l) . it is convenient to pick the gamma matrix basis as T'' = — CT3, F"^ = iai and F* = —a2- We 
can write explicitly: 



4(u. + ^) + £A, ^ac-™ il- 



= 0. (38) 



Note that the AdS2 radius R2 enters for the (r, () directions. At the AdS2 boundary, ( 0, the Dirac equation to 
the leading order is given by 




(39) 



The solution to this equation is given by the scaling function F*^') = Ae^Q + Be^Q'^ where e± are the real 
eigenvectors of U and the exponent is 




/ 

The conformal dimension of the operator O in the IR CFT is (5; = i + i^j . Comparing eq. (PO)) to the expression for 
the scaling exponent in Q , we conclude that the scaling properties and the AdS2 construction are unmodified by the 
magnetic field, except that the scaling exponents are now fixed by the Landau quantization. This " quantization rule" 
was already exploited in [22| to study de Haas- van Alphen oscillations. 



IV. SPECTRAL FUNCTIONS 



In this section we will explore some of the properties of the spectral function, in both plane wave and Landau 
level basis. We first consider some characteristic cases in the plane wave basis and make connection with the ARPES 
measurements . 
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A. Relating to the ARPES measurements 



In reality, ARPES measurements cannot be performed in magnetic fields so the holographic approach, allowing 
a direct insight into the propagator structure and the spectral function, is especially helpful. This follows from the 
observation that the spectral functions as measured in ARPES are always expressed in the plane wave basis of the 
photon, thus in a magnetic field, when the momentum is not a good quantum number anymore, it becomes impossible 
to perform the photoemission spectroscopy. 

In order to compute the spectral function, we have to choose a particular fermionic plane wave as a probe. Since 
the separation of variables is valid throughout the bulk, the basis transformation can be performed at every constant 
r-slice. This means that only the x and y coordinates have to be taken into account (the plane wave probe lives only 
at the CFT side of the duality). We take a plane wave propagating in the +x direction with spin up along the r-axis. 
In its rest frame such a particle can be described by 



^probc 



(41) 



Near the boundary (at oo) we can rescale our solutions of the Dirac equation making use of cqs. (|A23IIA24l IBll) : 



F, 



cP^(i) 



F 



\ 



-CP^(i) 



(42) 



with rescaled x defined after eq. (|A20[) . This representation is useful since we calculate the components £,±{ri,) related 
to the retarded Green's function in our numerics (we keep the notation of Q). 

Let Oi and Oi be the CFT operators dual to respectively Fi and Fi, and c^, Ck be the creation and annihilation 
operators for the plane wave state 4'probc- Since the states F and F form a complete set in the bulk, we can write 



4(0.) = 5^ ([/;,[/; 



(43) 



where the overlap coefl[icients Ui^uj) are given by the inner product between 4'probo and F: 



dxi^/ir°*pi.obc = - / dxe 



(44) 



with F = FHt'^, and similar expression for Ui involving ^_(rb). The constants Ui can be calculated analytically using 
the numerical value of ^±(rb), and by noting that the Hermite functions are eigenfunctions of the Fourier transform. 
We are interested in the retarded Green's function, defined as 



~i y d^dte'"*-'f^6l(t)Gg,(t,x) 

(0| [Oi{t,x),6i{0,0)] |0) 

Go 
Go 



(45) 



where Go is the retarded Green's function for the operator O. 

Exploiting the orthogonality of the spinors created by O and and using eq. (|43|). the Green's function in the 
plane wave basis can be written as 



Gf^(w,p,)=^tr 



U*,U*]G^ = 



(|C/i(Px)l'Gg,(w,0 + |?7i(p.)pGg,(w,0 



(46) 



In practice, we cannot perform the sum in eq. (|46p all the way to infinity, so we have to introduce a cutoff Landau 
level Zcut- In most cases we are able to make Icut large enough that the behavior of the spectral function is clear. 
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Using the above formalism, we have produced spectral functions for two different conformal dimensions and fixed 
chemical potential and magnetic field (Fig. 1). Using the plane wave basis allows us to directly detect the Landau 
levels. The unit used for plotting the spectra (here and later on in the paper) is the effective temperature Toff [H: 



(47) 



This unit interpolates between ^ at T//i = and T and is of or T/ fi — > oo. and is convenient for the reason that the 
relevant quantities (e. g., Fermi momentum) are of order unity for any value of fi and h. 






-0.6 -0.4 -0.2 0.0 0.2 0.4 




FIG. 1: Two examples of spectral functions in the plane wave basis for fi/T — 50 and h/T = 1. The conformal dimension is 
A = 5/4 (left) and A = 3/2 (right). Frequency is in the units of effective temperature Toff. The plane wave momentum is 
chosen to be A; = 1. Despite the convolution of many Landau levels, the presence of the discrete levels is obvious. 



B. Magnetic crossover and disappearance of the quasiparticles 



Theoretically, it is more convenient to consider the spectral functions in the Landau level basis. For definiteness 
let us pick a fixed conformal dimension A = | which corresponds to m = —j. In the limit of weak magnetic fields, 
h/T ^ 0, we should reproduce the results that were found in 

In Fig. (2. A) we indeed see that the spectral function, corresponding to a low value of /i/T, behaves as expected for 
a nearly conformal system. The spectral function is approximately symmetric about w = 0, it vanishes for |a;| < k, 
up to a small residual tail due to finite temperature, and for \uj\ ^ k it scales as o;^™. 

In Fig. (2.B), which corresponds to a high value of ^i/T, we see the emergence of a sharp quasiparticle peak. This 
peak becomes the sharpest when the Landau level I corresponding to an effective momentum kcs = y^2\qh\l coincides 
with the Fermi momentum kp. The peaks also broaden out when fcoff moves away from hp. A more complete view of 
the Landau quantization in the quasiparticle regime is given in Fig. IIVB( where we plot the dispersion relation {ui-k 
map). Both the sharp peaks and the Landau levels can be visually identified. 

Collectively, the spectra in Fig. 2 show that conformality is only broken by the chemical potential fj, and not by the 
magnetic field. Naively, the magnetic field introduces a new scale in the system. However, this scale is absent from 
the spectral functions, visually validating the the discussion in the previous section that the scale h can be removed 
by a rescaling of the temperature and chemical potential. 

One thus concludes that there is some value h'^ of the magnetic field, depending on /x/T, such that the spectral 
function loses its quasiparticle peaks and displays near-conformal behavior for h > h'^. The nature of the transition 
and the underlying mechanism depends on the parameters (/iq, T, A). One mechanism, obvious from the rescaling in 
eq. ([29| . is the reduction of the effective coupling q as h increases. This will make the influence of the scalar potential 
Aq negligible and push the system back toward conformality. Generically, the spectral function shows no sharp change 
but is more indicative of a crossover. 

A more interesting phenomenon is the disappearance of coherent quasiparticles at high effective chemical potentials. 
For the special case m = 0, we can go beyond numerics and study this transition analytically, combining the exact 
T = solution found in [2l| and the mapping ([50]) . In the next section, we will show that the transition is controlled 
by the change in the dispersion of the quasiparticle and corresponds to a sharp phase transition. Increasing the 
magnetic field leads to a decrease in phenomenological control parameter I'kF ■ This can give rise to a transition to a 
non- Fermi liquid when < 1/2, and finally to the conformal regime at h ~ h'^ when Vf^^ — and the Fermi surface 
vanishes. 
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FIG. 2: Some typical examples of spectral functions A{u!,kca) vs. oj in the Landau basis, fcoff = \/2\qh\n. The top 
four correspond to a conformal dimension A = | m = — i and the bottom four to A = | (m = 0). In each plot we 
show different Landau levels, labelled by index n, as a function of n/T and h/T. The ratios take values {fi/T, h/T) = 
(1, 1), (50, 1), (1, 50), (50, 50) from left to right. Conformal case can be identified when n/T is small regardless of h/T (plots in 
the left panel). Nearly conformal behavior is seen when both n/T and h/T are large. This confirms our analytic result that 
the behavior of the system is primarily governed by /i. Departure from the conformality and sharp quasiparticle peaks are seen 
when ji/T is large and h/T is small in (2.B) and (2.F). Multiple quasiparticle peaks arise whenever feeff = kp- This suggests 
the existence of a critical magnetic field, beyond which the quasiparticle description becomes invalid and the system exhibits 
a conformal- like behavior. As before, the frequency cu is in units of Tcff. 



12 




FIG. 3; Dispersion relation ui vs. kc« for /i/T = 50, h/T = 1 and A = | (m — —\)- The spectral function A{ijj,kc,s) is 
displayed as a density plot. (A) On a large energy and momentum scale, we clearly sees that the peaks disperse almost linearly 
(ai ~ Vpk), indicating that we are in the stable quasiparticle regime. (B) A zoom-in near the location of the Fermi surface 
shows clear Landau quantization. 




FIG. 4: Density of states D{uj) for m = — i and (A) ^/T = 50, h/T = 1, and (B) fi/T = 1, h/T = 1. Sharp quasiparticle peaks 
from the splitting of the Fermi surface are clearly visible in (A). The case (B) shows square-root level spacing characteristic of 
a (nearly) Lorentz invariant spectrum such as that of graphene. 

C. Density of states 

As argued at the beginning of this section, the spectral function can look quite different depending on the particular 
basis chosen. Though the spectral function is an attractive quantity to consider due to connection with ARPES 
experiments, we will also direct our attention to basis-independent and manifestly gauge invariant quantities. One of 
them is the density of states (DOS), defined by 

D{oj) = Y,A{uJ), (48) 
I 

where the usual integral over the momentum is replaced by a sum since only discrete values of the momentum are 
allowed. 

In Fig. 4, we plot the density of states for two systems. We clearly see the Landau splitting of the Fermi surface. 
A peculiar feature of these plots is that the DOS seems to grow for negative values of uj. This, however, is an artefact 
of our calculation. Each individual spectrum in the sum cq. (|48| has a finite tail that scales as w^™ for large uj, 
so each term has a finite contribution for large values of ui. When the full sum is performed, this fact implies that 
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lim D{uj) — >■ oo. The relevant information on the density of states can be obtained by regularizing the sum, which 

a;— f C30 

in practice is done by summing over a finite number of terms only, and then considering the peaks that lie on top of 
the resulting finite-sized envelope. The physical point in Fig. 4A is the linear spacing of Landau levels, corresponding 
to a non-relativistic system at finite density. This is to be contrasted with Fi g. 4 B where the level spacing behaves 
as oc -v/X, appropriate for a Lorentz invariant system and realized in graphene |24| . 



V. FERMI LEVEL STRUCTURE AT ZERO TEMPERATURE 

In this section, we solve the Dirac equation in the magnetic field for the special case to = (A = |). Although 
there are no additional symmetries in this case, it is possible to get an analytic solution. Using this solution, we 
obtain Fermi level parameters such as kp and vp and consider the process of filling the Landau levels as the magnetic 
field is varied. 



A. Dirac equation witli m = Q 

In the case to = 0, it is convenient to solve the Dirac equation including the spin connection (eq. IA2p rather than 
scaling it out: 



-^<y Or )^=(J^[UJ + qAt) + ! cr -UJift 



cr^^iOifx - cr^^i^yfy - Ai j ® 1 f ) = 0, (49) 



where A; = ^/2\qh\l are the energies of the Landau levels I ~ 0,1, ga = g^x = Qyy, ^tii") is given by eq. ([3]), and 
the gamma matrices are defined in eq. (jA12p . In the basis of eq. (|A12[) . the two components tpi and -02 decouple. 
Therefore, in what follows we solve for the first component only (we omit index 1). Substituting the spin connection, 
we have p^ : 



a'dr - ^(J'ii^ + qAt) -cT'^i3 + ^)-Xi]^ = 0, (50) 



i?2 V7 2i?2 ^ 2/ 

with tp = {yi,y2)- It is convenient to change to the basis 

m ) V -i 1 / V 2/2 



(51) 



which diagonalizes the system into a second order differential equation for each component. We introduce the dimen- 
sionless variables as in eqs. (|15m7p . and make a change of the dimensionless radial variable: 

(52) 



1-2' 



with the horizon now being at z = 0, and the conformal boundary at z = 1. Performing these transformations in eq. 
the second order differential equations for yi reads 



,2 , . 3/ , ,,,, , 15/ , 3/' , /" 



fdt + (— ^ + no.. 



l-z ' ' 4(1 - z)2 2(1 -z) 4 

+ i((w + (7^2)±^)2-zg;,-Af )^i=0, (53) 



The second component y^ obeys the same equation with [i i-> 
At T = 0, 



/ = 3z2(z-zo)(-2-^o), 20 = ^(4 + 1^2). (54) 
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The solution of this fermion system at zero magnetic field and zero temperature T = has been found in [2i| . To 
solve eq. (j53p , we use the mapping to a zero magnetic field system eq. (|29p . The combination = at non-zero h 
maps to /Xq^off = McffZcff at zero h as follows: 



/ ^2 / ^2 

fig^qJl- Q2^jj2 ■ 9fVQ^+H^ = VSqgpyl - — = /ig,off (55) 



where at T = we used + = 3. We solve eq. ([55|) for zero modes, i. c. uj ~ 0, and at the Fermi surface X — k, 
and implement eq. (j55p . 



Near the horizon (z = 0, / = 6z ), we have 



3 {P-q.cffY ,„2 



6z^yi'.2 + 12zy'i^2 + + - ^fJ Via = 0, (56) 

which gives the following behavior: 

yi;2-^"^^"N (57) 

with the scaling exponent v following from cq. (1321) : 

^=^V'6fc2-(M,,c//)2, (58) 

at the momentum k. Using Maple, we find the zero mode solution of eq. ([53]) with a regular behavior z~5+'' at the 
horizon [ll,[2l|: 

= A^,(.-l)i.-H^(z-z-o)--'^fi^) 

\^2 3 ^ 6 3zo(-z - -zo) y 

and 

^ i(-l+V2^,,o///zo) 

,-,(0) - 



y^*^^ = 7V2(^"l)^z-^+''(z-zo)-^-'^ 

X 2flU+^+4rM-?,eff.'^-^^^.l + 2i^: q / -x h 60 

y2 3 6 3zo(2 - zq) j 

where is the hypergeometric function and Ni,N2 are normalization factors. Since normalization factors are 
constants, we find their relative weight by substituting solutions given in eq. (|59p back into the first order differential 
equations at z ^ 0, 



A^2 ^ \/6fc \zo 



Mg,o///\/2zo 



(61) 



The same relations are obtained when calculations are done for any z. The second solution f]f^2-: with behavior z 2 
at the horizon, is obtained by replacing — > — in eq. (j59p . 

To get insight into the zero-mode solution (|59|) . we plot the radial profile for the density function ip^'^^'^ ip^'^^ for 
different magnetic fields in Fig. ([5]). The momentum chosen is the Fermi momentum of the first Fermi surface (see 
the next section). The curves are normalized to have the same maxima. Magnetic field is increased from right to left. 
At small magnetic field, the zero modes are supported away from the horizon, while at large magnetic field, the zero 
modes are supported near the horizon. This means that at large magnetic field the infiuence of the black hole to the 
Fermi level structure becomes more important. 
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FIG. 5: Density of the zero mode ip^^ip'^ vs. the radial coordinate z (the horizon is at 2; = and the boundary is at z = 1) for 
different values of the magnetic field h for the first (with the largest root for kp) Fermi surface. We set gp = 1 (h ^ H) and 

q = ^ (^,,cff 15y^l - -^1^). From right to left the values of the magnetic field are H = {0, 1.40, 1.50, 1.60, 1.63, 1.65, 1.68}. 
The amplitudes of the curves are normalized to unity. At weak magnetic fields, the wave function is supported away from the 
horizon while at strong fields it is supported near the horizon. 



B. Magnetic effects on the Fermi momentum and Fermi velocity at T = 

In the presence of a magnetic field there is only a true pole in the Green's function whenever the Landau level 
crosses the Fermi energy [22| 

2l\qh\ = fc|. (62) 

As shown in Fig. 2, whenever the equation is satisfied the spectral function A{u!) has a (sharp) peak. This is 
not surprising since quasiparticles can be easily excited from the Fermi surface. From eq. (j62p . the spectral function 
A{ij) and the density of states on the Fermi surface D{uj) are periodic in with the period 

where Ap = nkj^ is the area of the Fermi surface (22j . This is a manifestation of the de Haas- van Alphen quantum 
oscillations. At T = 0, the electronic properties of metals depend on the density of states on the Fermi surface. 
Therefore, an oscillatory behavior as a function of magnetic field should appear in any quantity that depends on the 
density of states on the Fermi energy. Magnetic susceptibility [22j and magnetization together with the supercon- 
ducting gap [2^ have been shown to exhibit quantum oscillations. Every Landau level contributes an oscillating term 
and the period of the Z-th level oscillation is determined by the value of the magnetic field h that satisfies eq. (|62|) for 
the given value of fc^. Quantum oscillations (and the quantum Hall effect which we consider later in the paper) are 
examples of phenomena in which Landau level physics reveals the presence of the Fermi surface. The superconducting 
gap found in the quark matter in magnetic fields [2^ is another evidence for the existence of the (highly degenerate) 
Fermi surface and the corresponding Fermi momentum. 

Generally, a Fermi surface controls the occupation of energy levels in the system: the energy levels below the Fermi 
surface are filled and those above are empty (or non-existent). Here, however, the association to the Fermi momentum 
can be obscured by the fact that the fermions form highly degenerate Landau levels. Thus, in two dimensions, in 
the presence of the magnetic field the corresponding effective Fermi surface is given by a single point in the phase 
space, that is determined by np, the Landau index of the highest occupied level, i.e., the highest Landau level below 
the chemical potential ^. Increasing the magnetic field. Landau levels 'move up' in the phase space leaving only the 



^ We would like to thank Igor Shovkovy for clarifying the issue with the Fermi momentum in the presence of the magnetic field. 
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lower levels occupied, so that the effective Fermi momentum scales roughly (excluding interactions) as a square root 
of the magnetic field, kp ~ ^/np ~ k^""^ ^1 — h/hmax- High magnetic fields drive the effective density of the charge 
carriers down, approaching the limit when the Fermi momentum coincides with the lowest Landau level. 

Many phenomena observed in the paper can thus be qualitatively explained by Landau quantization. As discussed 
before, the notion of the Fermi momentum is lost at very high magnetic fields. In what follows, the quantitative Fermi 
level structure at zero temperature, described by fc^ and vp values, is obtained as a function of the magnetic field 
using the solution of the Dirac equation given by cqs. (j59l60|) . As in [ll|, we neglect first the discrete nature of the 
Fermi momentum and velocity in order to obtain general understanding. Upon taking the quantization into account, 
the smooth curves become combinations of step functions following the same trend as the smooth curves (without 
quantization). While usually the grand canonical ensemble is used, where the fixed chemical potential controls the 
occupation of the Landau levels j26| . in our setup, the Fermi momentum is allowed to change as the magnetic field is 
varied, while we keep track of the IR conformal dimension v. 



keff 




H 



0.5 1.0 1.5 




FIG. 6; Effective momentum fee// vs. the magnetic field h ^ H (we set gp = I, q = -^). As we increase magnetic field the 
Fermi surface shrinks. Smooth solid curves represent situation as if momentum is a continuous parameter (for convenience), 
stepwise solid functions are the real Fermi momenta which are discretized due to the Landau level quantization: kp — > \/2\qh\l 
with I = 1,2, . . . where ^J2\qh\l axe Landau levels given by dotted lines (only positive discrete kp are shown). At a given h 
there are multiple Fermi surfaces. From right to left are the first, second etc. Fermi surfaces. The dashed-dotted line is vi^^ — 
where kp is terminated. Positive and negative keff correspond to Fermi surfaces in two components of the Green's function. 



The Fermi momentum is defined by the matching between IR and UV physics Q, therefore it is enough to know 
the solution at cj = 0, where the matching is performed. To obtain the Fermi momentum, we require that the zero 
mode solution is regular at the horizon (tp^^") ^ z^^^^) and normalizable at the boundary. At the boundary z 1, 
the wave function behaves as 



1 \ , , /, / 



a(l-z)^-'"(^Q j +&(l-z)5+"(^:j' j. (64) 
To require it to be normalizable is to set the first term a = 0; the wave function at 2 ~ 1 is then 

^(0) ^ (l-z)i+™ (^5^ . (65) 

Eq. (|55|) leads to the condition lim2_j.i(z — l)^^^'^{yf^ + i^i^') = 0, which, together with eq. ([5^ . gives the following 
equation for the Fermi momentum as function of the magnetic field [isl . [2]| 

2Fi(l + 1 + 21., 1(1 - *y2)) _ %i^-iti,,.ff ^gg^ 



2Fi(i/ + '-i^, \ + 1 + 2v, 1(1 - i^/2)) kp{-2i + ^/2) ' 

with V = given by eq. (j58p . Using Mathematica to evaluate the hypergeometric functions, we numerically solve 
the equation for the Fermi surface, which gives effective momentum as if it were continuous, i.e. when quantization 
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FIG. 7: Landau level numbers n coresponding to the quantized Fermi momenta vs. the magnetic field h ^ H for the three 
Fermi surfaces with positive kp. We set gp ~ 1, q = From right to left are the first, second and third Fermi surfaces. 



is neglected. The solutions of eq. (|66)) are given in Fig. ([6]). There are multiple Fermi surfaces for a given magnetic 
field h. Here and in all other plots we choose gp = 1, therefore h ^ H, and q = In Fig.(l6|), positive and negative 
kp correspond to the Fermi surfaces in the Green's functions Gi and G2. The relation between two components is 
G2{(jJ,k) = Gi{uj, —k) 0, therefore Fig.® is not symmetric with respect to the x-axis. Effective momenta terminate 
at the dashed line Vkp = 0. Taking into account Landau quantization of kp — ^ ^J2\qh\l with I = 1,2..., the plot 
consists of stepwise functions tracing the existing curves (we depict only positive kp). Indeed Landau quantiaztion 
can be also seen from the dispersion relation at Fig. (3), where only discrete values of effective momentum are allowed 
and the Fermi surface has been chopped up as a result of it Fig. (3B). 

Our findings agree with the results for the (largest) Fermi momentum in a three-dimensional magnetic system 
considered in [23l, compare the stepwise dependence kp{h) with Fig. (5) in j27j . 

In Fig. ([7]), the Landau level index / is obtained from kpih) = ^J2\qh\l where kp{h) is a numerical solution of eq. 
Only those Landau levels which are below the Fermi surface are filled. In Fig.®, as we decrease magnetic field 
first nothing happens until the next Landau level crosses the Fermi surface which corresponds to a jump up to the 
next step. Therefore, at strong magnetic fields, fewer states contribute to transport properties and the lowest Landau 
level becomes more important (see the next section). At weak magnetic fields, the sum over many Landau levels has 
to be taken, ending with the continuous limit as /i — > 0, when quantization can be ignored. 

In Fig. ([8|), we show the IR conformal dimension as a function of the magnetic field. We have used the numerical 
solution for kp. Fermi liquid regime takes place at magnetic fields h < he, while non- Fermi liquids exist in a narrow 
band at he < h < h'^, and at h'^ the system becomes near-conformal. 

In this figure we observe the pathway of the possible phase transition exhibited by the Fermi surface (ignoring 
Landau quantization): it can vanish at the line vj^^ = 0, undergoing a crossover to the conformal regime, or cross 
the line v^f =1/2 and go through a non- Fermi liquid regime, and subsequently cross to the conformal phase. Note 
that the primary Fermi surface with the highest kp and Vkp seems to directly cross over to conformality, while the 
other Fermi surfaces first exhibit a " strange metal" phase transition. Therefore, all the Fermi momenta with i>kp > 
contribute to the transport coefficients of the theory. In particular, at high magnetic fields when for the first (largest) 

Fermi surface k^p^ is nonzero but small, the lowest Landau level n = becomes increasingly important contributing 
to the transport with half degeneracy factor as compared to the higher Landau levels. 

In Fig. ini we plot the Fermi momentum kp as a function of the magnetic field for the first Fermi surface (the 
largest root of eq. ^^). Quantization is neglected here. At the left panel, the relatively small region between the 
dashed lines corresponds to non- Fermi liquids < ly < ^. At large magnetic field, the physics of the Fermi surface 
is captured by the near horizon region (see also Fig. ([5])) which is AdS2 x M^. At the maximum magnetic field, 
Hmax = \/3 ~ 1.73, when the black hole becomes pure magnetically charged, the Fermi momentum vanishes when 
it crosses the line i^kp = 0. This only happens for the first Fermi surface. For the higher Fermi surfaces the Fermi 
momenta terminate at the line = 0, Fig. Note the Fermi momentum for the first Fermi surface can be almost 

fully described by a function kp — fc^'^^wl — It is tempting to view the behavior kp ^ \/Hmax — H as a phase 
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FIG. 8: Left panel. The IR conformal dimension v = i/^p calculated at the Fermi momentum vs. the magnetic field h ^ H 
(we set (/F = l, q = ■^)- Calculations are done for the first Fermi surface. Dashed line is for = \ (at He = 1.70), which is the 

border between the Fermi liquids v > ^ and non- Fermi liquids v < ^. Right panel. Phase diagram in terms of the chemical 
potential and the magnetic field /i^ + /i^ = 3 (in dimensionless variables h = QfH, /i = QfQ', we set qf = 1). Fermi liquids are 
above the dashed line (H < He) and non-Fermi liquids are below the dashed line {H > He). 



transition in the system although it strictly follows from the linear scaling for _ff = by using the mapping (j29p . 
(Note that also fi ~ gpQ = pfVS — H"^.) Taking into account the discretization of /cf, the plot will consist of an 
array of step functions tracing the existing curve. Our findings agree with the results for the Fermi momentum in a 
three dimensional magnetic system considered in (27j . compare with Fig. (5) there. 
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FIG. 9: Fermi momentum kp vs. the magnetic field h ^ H (we set gp = 1, q ~ -^) for the first Fermi surface. Left panel. The 
inner (closer to x-axis) dashed line is Vkp = and the outer dashed line is Ukp = \ , the region between these lines corresponds 
to non-Fermi liquids < vtp < |. The dashed-dotted line is for the first Landau level fci — y/2qH. The first Fermi surface 
hits the border-line between a Fermi and non-Fermi liquids v = ^ at He ~ 1.70, and it vanishes at Hmax = y/S = 1.73. Right 

panel. Circles are the data points for the Fermi momentum calculated analytically, solid line is a fit function fc™"^ \J ^ ~ ^ 
with fc^"^ = 12.96. 



The Fermi velocity given in eq. ((27)) is defined by the UV physics; therefore solutions at non-zero uj are required. 
The Fermi velocity is extracted from matching two solutions in the inner and outer regions at the horizon. The Fermi 
velocity as function of the magnetic field for > ^ is [3, Hlj 



..(0) 

= hm ^1 , + , (67) 

where the zero mode wavefunction is taken at kp ea.(|59p. 

We plot the Fermi velocity for several Fermi surfaces in Fig. [TU] and for the first Fermi surface in Fig. [TT] 
Quantization is neglected here. The Fermi velocity is shown for v > ^. It is interesting that the Fermi velocity 
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vanishes when the IR conformal dimension is i^kp = \- Formally, it follows from the fact that vp ^ {2v — 1) [7|. 
The first Fermi surface is at the far right. Positive and negative vp correspond to the Fermi surfaces in the Green's 
functions Gi and G2, respectively. The Fermi velocity vp has the same sign as the Fermi momentum kp. At small 
magnetic field values, the Fermi velocity is very weakly dependent on H and it is close to the speed of light; at large 
magnetic field values, the Fermi velocity rapidly decreases and vanishes (at He = 1.70 for the first Fermi surface 
Fig. ()lip ). Geometrically, this means that with increasing magnetic field the zero mode wavefunction is supported 
near the black hole horizon fig.®, where the gravitational redshift reduces the local speed of light as compared to 
the boundary value. It was also observed in [tI l2l| at small fermion charge values. 
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FIG. 10: Fermi velocity vf vs. the magnetic field h ^ H (we set gp = 1, q = -^) for the regime of Fermi liquids v > \. Fermi 

velocity vanishes at Ukp = | (x-axis). The multiple lines are for various Fermi surfaces in ascending order, with the first Fermi 
surface on the right. The Fermi velocity vf has the same sign as the Fermi momentum kp- As above, positive and negative 
Vp correspond to Fermi surfaces in the two components of the Green's function. 



vF 

1.0 



0.8 
0.6 
0.4 
0.2 



""O.O 0.5 1.0 15 2.0" 

FIG. 11: Fermi velocity vp vs. the magnetic field h ^ H (we set gp = 1, q = -^) for the first Fermi surface. Fermi velocity 
vanishes at = i at He ~ 1.70. The region H < He corresponds to the Fermi liquids and quasiparticle description. 



VI. HALL AND LONGITUDINAL CONDUCTIVITIES 

In this section, we calculate the contributions to Hall a^y and the longitudinal a^x conductivities directly in 
the boundary theory. This should be contrasted with the standard holographic approach, where calculations are 
performed in the (bulk) gravity theory and then translated to the boundary field theory using the AdS / CFT dictionary. 
Specifically, the conductivity tensor has been obtained in [l^ by calculating the on-shell renormalized action for the 
gauge field on the gravity side and using the gauge/gravity duality Am ^ j/i to extract the R charge current-current 
correlator at the boundary. Here, the Kubo formula involving the current-current correlator is used directly by 
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utilizing the fermion Green's functions extracted from holography in [7|. Therefore, the conductivity is obtained for 
the charge carriers described by the fermionic operators of the boundary field theory. 

The use of the conventional Kubo formulo to extract the contribution to the transport due to fermions is validated 
in that it also follows from a direct AdS/CFT computation of the one-loop correction to the on-shell renormalized 
AdS action [Tst - We study in particular stable quasiparticles with v > ^ and at zero temperature. This regime 
effectively reduces to the clean limit where the imaginary part of the self-energy vanishes ImS — > 0. We use the 
gravity- "dressed" fermion propagator from eq. (|27p and to make the calculations complete, the "dressed" vertex is 
necessary, to satisfy the Ward identities. As was argued in p^ . the boundary vertex which is obtained from the bulk 
calculations can be approximated by a constant in the low temperature limit. Also, according to [25|, the vertex 
only contains singularities of the product of the Green's functions. Therefore, dressing the vertex will not change 
the dependence of the DG conductivity on the magnetic field [l^ . In addition, the zero magnetic field limit of the 
formulae for conductivity obtained from holography [l5| and from direct boundary calculations [l8| are identical. 



A. Integer quantum Hall effect 

Let us start from the "dressed" retarded and advanced fermion propagators 0: Gr is given by eq. (|27p and 
Ga = G^. To perform the Matsubara summation we use the spectral representation 

A [ duj A(uj,k) 

G{tLJ„,k) = / 68 

J It: u — iLOn 

with the spectral function defined as A{uj, k) = — ^ImG ft juj , k) = ^{Gji{oj, k) — GA{^, k)). Generalizing to a non-zero 
magnetic field and spinor case [1^, the spectral function [2^ is 

Aiu., k) = ie-^ ti-lYi-h^v.) ( . w , + ^ -i^o) , (69) 

TT ^ + ep + i-'i{i^,kp) - EiY + Y.2[i^,kFY j 



where sp — vpkp is the Fermi energy, Ei ~ vp^2\qh\l is the energy of the Landau level, j(k) = P^Li{j^) — 

P-i_L/_i(-^^) with spin projection operators P± = (l±i7^7^)/2, we take c = 1, the generalized Laguerre polynomials 

are L^(z) and by definition £„(z) — L'^{z), (we omit the vector part fc-f, it does not contribute to the DG conductivity), 
all 7's are the standard Dirac matrices, vp and kp are real constants (we keep the same notations for the constants 
as in 0). The self-energy E ~ uY'^''f contains the real and imaginary parts, S = Si + iY,2- The imaginary part 
comes from scattering processes of a fermion in the bulk, e. g. from pair creation, and from the scattering into the 
black hole. It is exactly due to inelastic/dissipative processes that we are able to obtain finite values for the transport 
coefficients, otherwise they are formally infinite. 

Using the Kubo formula, the DG electrical conductivity tensor is 

ImH,^ 

a^AVt) = lim ^, (70) 

where Hijii^m O + iO+) is the retarded current-current correlation function; schematically the current density 
operator is j^(t, x) ~ qvp i^aiT, x)^''"ipa{T, x). Neglecting the vertex correction, it is given by 

[lilrn) = q^vlT j T^tr (7'G(zc^„, fc)7^G(^a;„ + lU^, k)) . (71) 

71 — — oo ' 

The sum over the Matsubara frequency is 



1 1 niixix) — n(uj2) 

Taking iflm — > -I- iO^ , the polarization operator is now 



jnV^ ^ ^ ^ n[0Ji)-n[uj2) ^^^^ 



TT (r>\ c?wi da;2 nFD(wi) - nFD(w2) f d'^k . , 

n + .,-u.2 J (^tr(7^(-i,^')y^(-2,fc)), (73) 
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where the spectral function A^lo, k) is given by eq. (j69p and n-pD{uj) is the Fermi-Dirac distribution function. Evalu- 
ating the traces, we have 

a., = ^ ^^'^FihiVF^qhl ^^ J. {s,,{S,k-i + ^^-1,.) + .6.,,-sgn(gfe)(J,fc_i - S,^,,,)} 

l,k=0 

— - (tanh — ^ — tanh — ^ ) 
27r ^ 2T 2T' 



(cDi-£;,)2+SlK) ' ^ ' 'V V('^2-i?fc)2 + S2(c.2) 

with = i^i + We have also introduced uji-2 = '^i;2 + ep + Si(a;i;2) with Sij being the antisymmetric tensor 
(ei2 = 1), and I]i;2(ct') = Si:2(a;, kp). In the momentum integral, we use the orthogonality condition for the Laguerre 
polynomials L dxe'^ Li{x)L}^{x) = (5;^. 

From eq. (|74p . the term symmetric/ antisymmetric with respect to exchange cji -(->■ W2 contributes to the diagonal/off- 
dialgonal component of the conductivity (note the antisymmetric term tifdIi^i) ~ '^fd('^2))- The longitudinal and 
Hall DC conductivities (O — > 0) are thus 

2q^{hiVFf\qh\ f°° dio 



ttT 7-00 2^ cosh^ ^ 



X 



E ( 7^ ^ + ^ ) ( 7^ F !^^v2^ ^ + (^'+1 ^ ) ' (75) 

^ \ [uj - El)- + T,^{uj) J \{uj - Ei+iY + T.2{uj) J 



1=0 

q^{hiVFfsgn{qh) 
i^h, 



Vh 



/OO ^ ^ / 1 \ 



where a) = w + e^? + I]i(aj)). The filling factor Vh is proportional to the density of carriers: = |^^|^^^^ n (we derive 
this relation below eq. (HU). The degeneracy factor of the Landau levels is a/: ao = 1 for the lowest Landau level 
and Of; = 2 for Z = 1, 2 . . . . Substituting the filling factor Vh back to eq. ([76|) . the Hall conductivity can be written as 

f^xy = j^, (77) 



where p is the charge density in the boundary theory, and both the charge q and the magnetic field h carry a sign (the 
prefactor {—hivp) comes from the normalization choice in the fermion propagator eqs. (j27l69p as given in 0, which 
can be regarded as a factor contributing to the effective charge and is not important for further considerations) . The 
Hall conductivity eq. ([77)) has been obtained using the AdS / CFT duality for the Lorentz invariant 2 + 1-dimensional 
boundary field theories in [l^ . We recover this formula because in our case the translational invariance is maintained 
in the x and y directions of the boundary theory. 

Low frequencies give the main contribution in the integrand of eq. (|76p . Since the self-energy satisfies Si(a;) 
T,2{ijj) ^ iiJ^" and we consider the regime v > |; we have Si S2 ^ at w ^ (self-energy goes to zero faster 
than the to term). Therefore, only the simple poles in the upper half-plane wq = —£f ± £-/ + Si -|- 1^2 contribute to 
the conductivity where Ei ^ E2 ~ (— Ef ± Eif''' are small. The same logic of calculation has been used in [2^. We 
obtain for the longitudinal and Hall conductivities 

2q^{hiVF)^^2 ( 1 1 + cosh ^ cosh f- ' 



X 



ttT \^l + cosh^ ^ (cosh^ -hcoshf-)2 

q^{hiVFffigTL{qh) 



'xy 



2 (tanh g + f (tanh + tanh £^)) , 



(79) 



where the Fermi energy is £f = vpkp and the energy of the Landau level is Ei = VFy'2\qh\l. Similar expressions were 
obtained in [2^. However, in our case the filling of the Landau levels is controlled by the magnetic field h through 
the field-dependent Fermi energy VF{h)kF{h) instead of the chemical potential fi. 



22 



At T = 0, cosh y — )• and tanh^ = 1 — 2nFD('^) ^ sgnw. Therefore the longitudinal and Hall conductivities 



are 



■2^«0eF,Bi= ^; Xn-OeF,-E„, (80) 



(=1 



q^ihivp) sgii{qh) 



n 

e{h,v,fsg.[qK^ X 2(1 + 2n)^(s. - E^)eiE^,, ^ e.), (81) 



2 l + 2^0(ej.-SO 



1=1 



where the Landau level index runs n — 0,1^ ... . It can be estimated as n 



when vp 7^ ([ ] denotes the 



_2|9/i|_ 

integer part), with the average spacing between the Landau levels given by the Landau energy Ui?-\/2|g/i|. Note that 
ep = epih). We can see that eq. (|8T1) expresses the integer quantum Hall effect (IQHE). At zero temperature, as we 
dial the magnetic field, the Hall conductivity jumps from one quantized level to another, forming plateaus given by 
the filling factor 

J.,, = ±2(l + 2n)=±4(?7,+ i), (82) 

with n = 0,1,.... (Compare to the conventional Hall quantization Vh = ±4n, that appears in thick graphene). 
Plateaus of the Hall conductivity at T = follow from the stepwise behavior of the charge density p in cq. (|77|) : 

p ^ 4(n + ^)e{eF - ErMEn+i - sp), (83) 

where n Landau levels are filled and contribute to p. The longitudinal conductivity vanishes except precisely at the 
transition point between the plateaus. In Fig. I12[ we plot the longitudinal and Hall conductivities at T = 0, using only 
the terms after x sign in eq. ((79| . In the Hall conductivity, plateau transition occurs when the Fermi level (in Fig. 
[T2|) of the first Fermi surface ep = vp(h)kp{h) (Figs. I9|lip crosses the Landau level energy as we vary the magnetic 
field. By decreasing the magnetic field, the plateaus become shorter and increasingly more Landau levels contribute 
to the Hall conductivity. This happens because of two factors: the Fermi level moves up and the spacing between the 
Landau levels becomes smaller. This picture does not depend on the Fermi velocity as long as it is nonzero. 

In the boundary field theory, wc express the charge density of the carriers (difference between the densities of 
"electrons" and "holes") through the Fermi energy £p: 

n = tr f7°G(r, 0)) , r ^ 0, (84) 



where G{t,x) is the translation- invariant part of the Green's function G{t,x) from eq. (|68p f26| . Using the spectral 
function representation eq. (|69p . the charge density reads 



We express the Matsubara sum in terms of the contour integral over real frequencies: 

T f2 F{ioj^)~^-^ I d^tanh^F(z), (86) 

where C runs ant i- clockwise and encircles the poles of tanh along the upper and lower half imaginary axis. We have 
for the charge density 

- = 2 y (^ y_ ^ -tr (/A(., k)) . (87) 
Substituting the spectral function eq. (|69p and integrating over momenta, we obtain 

2\qh\h'iVp f°° dui , uj „ , ^ /I „ 
^ = - ' ' / — tanh— EsHE^MT^ + 
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where the degeneracy factor is ao = 1 for the lowest Landau level and ai = 2 for the higher Landau levels I > 1, 
uj = uj + £p + Ei(u;). Integrating over frequencies and taking into account that E2 is effectively very small near the 
Fermi surface, we obtain 



X 2 tanh + 5:(tanh '-^^ + tanh '-^) . (89) 



2T 2T 

\ 1=1 

Comparing with eq. (|79p . we obtain the relation = \qh\hiVF When the Fermi energy vanishes {ep = 0), the 
spectral function eq. ()69p is even in w. From eq. (|88p . the carrier density of stable quasiparticles vanishes when 
£p = 0. At the end of this section, we discuss a situation with no stable charge carriers and physical consequences of 
it. 




FIG. 12: Hall conductivity a^y and longitudinal conductivity axx vs. the magnetic field h ~> H a,t T — (we set gp = 1, 
1 ~ -^)- Contribution from the first Fermi surface is taken. By decreasing the magnetic field, the Fermi surface crosses the 
Landau levels producing the Hall conductivity plateaus characteristic for IQHE. Longitudinal conductivity has picks at the 
beginning of each plateau. The right panel is a zoom-in for large h of the left one. 



Eqs. ()79ll89p are obtained assuming that the states are localized around the Landau levels. In QHE models, 
impurities are added to prevent the states from 'spilling' between the Landau levels and to provide the necessary 
occupation number of the levels. In our holographic calculations, however, the complex self-energy arises not from the 
impurities but from various scattering processes into the black hole. Here, the limit ImE — )■ has been considered, 
which corresponds to a simplified field theory model [2(| (the cited reference also considers the case with impurities) . 
This approximation suffices to obtain the integer QHE [2g and for our initial studies of the fractional QHE. We leave 
the implementation of a physical model with impurities for future work. 



B. Fractional quantum Hall effect 

In a holographic setting, using the AdS geometry is equivalent to a calculation in a box. Therefore, for large 
enough fermion charge q, there are multiple Fermi surfaces as shown in Figs. I6I10I Labelling the Fermi surfaces with 
> ^ by m = 1, 2, . . . , we represent the spectral function A{uj, k) as a sum over the spectral functions of individual 



Fermi surfaces given by eq. ([M)) [2l| . Ignoring the mixing term, the DC conductivity becomes a direct sum over the 
individual conductivities. By decreasing the magnetic field, new Fermi surfaces gradually appear as can be seen in 
Figs. I6I7I Therefore the conductivity tensor is 

^.,=Y.^t^(^^h^^a.-h). (90) 

m 

where el™'' involves the Fermi momentum /c^"* and velocity w^'' , respectively; at the maximum magnetic field h^max 
a new fc^-* opens up; h'^max is found numerically. 
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Including one, two, three and four Fermi surfaces, we obtain the following quantization rule for the filling factor in 
the Hall conductivity: 

1 FS : t^/i = 2(1 + 27i), plateaus -> 2, 6, 10, . . . , 
2FS's: i^/i = 4(1 + n + fc), plateaus ^ 4, 8, 12, . . . , 
3FS's: i^h^2{3 + 2{n + k + p)), plateaus ^ 6, 10, 14, ... , 

4FS's: i^/i = 4(2 + 77, + fc+p + r), plateaus 8, 12, 16, ... , (91) 

with 77, r = 0, 1, . . . . An odd number of Fermi surfaces produces the plateaus present in the IQHE, while an even 
number of Fermi surfaces produces the additional plateaus appearing in the fractional QHE (FQHE). For a large 
enough fermion charge q, many Fermi surfaces contribute, and the primary effect of the change in H is the opening 
of a new Fermi surface, rather than the occupation of the next plateau. Thus at large q we expect a filling fraction 
pattern at large h to become 

i^h = ±2j, (92) 

where j = 1,2, . . . is the effective Landau level index counting the number of contributing Landau levels. This is 
indeed observed in the FQHE at strong magnetic fields. The quantization rule (|91]) persists as long as new Fermi 
surfaces open up with decreasing h. However, the first two plateaus present in the FQHE ~ 0,±1 are absent in 
eq. ([M)) . In order to get the Hall plateau i^h ~ ±1, the mixing term between two Fermi surfaces should probably be 
taken into account (incoherent superposition), whereas the conductivity (|90l) includes the diagonal terms only. We 
discuss the issue with Vh =Q further. 

In Fig. we plot the Hall and longitudinal conductivities at T = with three Fermi surfaces contributing (eq. 

IM]) . where the individual conductivities cr^™-' are given by eq. ([75]). We fit the Fermi momenta by 



4™^ = 4™L\/l-y + ^^™\ (93) 



with = 12-96, <S(i' = 0., = 10-29, ^(2) = 1.5, kf^^^ = 9.75, 5^^^ = 3, and use eq. ^ together 

with the numerical solutions for v'^^^ in Fig. [131 In Fig. [131 at strong magnetic fields, the Hall conductivity plateau 
Vh ~ 4 originates from two Fermi surfaces together with the plateaus Vh ~ and Vh = 6 when one and three Fermi 
surfaces contribute, respectively. As we decrease the magnetic field further, three Fermi surfaces produce plateaus 
characteristic for IQHE eq. The longitudinal conductivity shows a Dirac delta-like peak at the beginning of 

each plateau. Since a finite contribution to the conductivity arises as one of the three Fermi surfaces crosses the 
next Landau level, the pattern is less regular (i. e., the plateaus have changing length) than in the case when only 
one Fermi surface contributes. In Fig. I13[ we compare the Hall conductivities with one and three Fermi surfaces 
participating. The irregular behavior of the Hall conductivity is explained naturally from the picture with multiple 
Fermi surfaces. Qualitatively similar regularity of the plateaus' length is seen in experiments on thin films of graphite 
at strong magnetic fields The actual physics behind, however, might be quite different as in this system multiple 
sheets of the Fermi surface arise due to the (hexagonal) lattice on the UV scale, an effect which is beyond the scope 
of our current model. 

The somewhat regular pattern behind the irregular behavior can be understood as a consequence of the appearance 
of a new energy scale: the average distance between the Fermi levels. For the case of Fig. [T31 we estimate it to be 
< e^'' — e^P'^^'^ >= 4.9 with ttt, = 1, 2. The authors of [2^ explain the FQHE through the opening of a gap in the 
quasiparticle spectrum, which acts as an order parameter related to the particle-hole pairing and is enhanced by the 
magnetic field (magnetic catalysis). Here, the energy gap arises due to the participation of multiple Fermi surfaces. 

A pattern for the Hall conductivity that is strikingly similar to Fig. ([13]) arises in the AA and AB-stacked bilayer 
graphene, which has different transport properties from the monolayer graphene [s^l, compare with Figs. (2,5) there. 
It is remarkable that the bilayer graphene also exhibits the insulating behavior in a certain parameter regime. This 
agrees with our findings of metal-insulating transition in our system. 



C. Metal-" strange metal" phase transition 



The previous discussion of conductivities and QHE is valid provided that the Fermi velocity is nonzero. However, 
we have shown that vp vanishes at relatively strong magnetic fields (for the first Fermi surface it happens at he as 
in Fig. [8l and [TT|) . In the AdS/CFT setting, the Fermi velocity vanishes when the IR anomalous dimension is = i 
signalling the onset of a nontrivial power law dispersion in Green's function G~^{u!) ^ to — Vfkj_ -^uP'^ (the pole in the 
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FIG. 13: Hall conductivity a^y and longitudinal conductivity cjxx vs. the magnetic field h H at T — (we set gp = 1, 
q = -^). Contribution from the first three Fermi surfaces are taken. At strong magnetic fields, the Hall conductivity plateau 
~ 4: appears from two Fermi surfaces together with plateaus Uh = 2 and 1/^ = 6 when one and three Fermi surfaces 
contribute, respectively. This quantization rule is characteristic for the FQHE. At intermediate and weak magnetic fields, 
the Hall conductivity plateaus are produced as one of the three Fermi surfaces crosses the Landau levels, resulting in the 
quantization rule of the IQHE. Irregular pattern in the length of the plateaus is observed in the experiment on thin films of 
graphite at strong magnetic fields [2j]. The right panel is a zoom-in for large h of the left one. 



self-energy E ^ oj^" and the pole in the prefactor of the linear term ~ w Q)- Vanishing of was observed 

in 21| at large enough fermion charge. Note that if is zero for some interval of the magnetic field, it leads to the 
Hall plateau with the filling factor i^h — present in FQHE. 

The vanishing of the Fermi velocity of the stable quasiparticle leads to zero carrier density at leading order: 

i)F = ^ n = 0. (94) 

This means that all contribution to conductivity comes from the other terms, containing the contribution from the 
non-Fermi liquid excitations and the conformal regime. This qualitatively changes the transport properties of the 
system, as can be seen in Fig. 1141 

The finite offset magnetic field has been observed in experiments on highly oriented pyrolitic graphite in magnetic 
fields [29| . In particular, analyzing the basal-plane resistivity an approximate scaling relation between the critical 
temperature of the metal-semiconducting transition and the magnetic field has been found Tc ^ \/h — he- It suggests 
that at T = 0, there is a threshold magnetic field he above which the resistivity qualitatively changes. Interestingly, 
the existence of such a threshold magnetic field follows from the AdS/CFT calculations {he when i/ = i). 

A phase transition is usually governed by an order parameter which exhibits a critical behavior. In our case, there is 
no such order parameter. However, it is interesting to note that the Fermi momentum, according to eq. (|93|) . behaves 
as fci^ '--^ y/hmax — h, which is in line with the postulated critical behavior in the system, while the Fermi surface itself 
behaves as order parameter. 

To obtain a complete picture of the metal-" strange metal" phase transition, one needs to perform calculations in 
the non-Fermi liquid regime taking into account the unstable quasiparticle pole. It is also necessary to study the 
temperature dependence of the DC conductivities aocyiT) and a^xiT). We leave it for the future study. 



VII. ABSENCE OF THE SIGN PROBLEM IN HOLOGRAPHY. 



In this section we show that the fermion determinant in the gravity dual theory docs not have a sign problem and 
hence can be simulated by a lattice Monte-Carlo algorithm. Up to recently, most of the work on AdS/CFT and applied 
holography focused on the classical gravity (leading in field theory) limit. However, many thermodynamic and 
electric properties depend on matter fields (e. g., the electrical conductivity depends on whether or not the theory has a 
Fermi surface). In classical gravity, the Einstein-Maxwell sector decouples, and matter fields run in loops representing 
quantum oscillations. In order to include matter fields in the bulk, one needs to calculate loop corrections, which 
corresponds to going beyond the leading order in 1/A'^. A study of one- loop bulk physics was done in |3l| and recently 
in p^ . It shows that analytical calculations of quantum corrections in the bulk are quite involved. The study of 
quantum oscillations in the gravity dual will likely improve our understanding of finite density systems in general. 
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FIG. 14: Comparison of the Hall conductivities a^y vs. the magnetic field h ^ H from one Fermi surface (dashed line) and 
from three Fermi surfaces (solid line). We set gp = 1 and q = At strong magnetic fields, a new plateau Vh = 4: 

appears in the multiple Fermi surface picture, yielding a pattern characteristic of FQHE. 



As is well known, a finite density field tlieory in most cases cannot be simulated on tlic lattice because of the 
infamous sign problem [l[. In the field theory action, chemical potential is introduced via the term '0/i7°'0 which 
is Hermitian and therefore gives a complex determinant. At the same time, in the bulk action, finite density is 
introduced through the electrically charge black hole, and docs not involve even matter fields. This is the reason 
why the applied holography gives universal predictions. In the leading order, the minimal gravitational dual at finite 
density and temperature is the electrically charged AdS-Rcissncr-Nordstrom black hole, where only the metric and 
Maxwell fields are present. Therefore, the Einstein-Maxwell sector gives results which do not depend, for example, on 
the charge and mass of matter fields in the gravitational bulk spacctime, i. e. are universal for a class of field theories 
with different charge and scaling dimensions of the operators. The fact that the chemical potential enters via the 
electric field in the covariant derivative leads to the real and positive definite fermion determinant which is suitable 
for lattice simulations. Wc show it formally below. 

In a semiclassical approach to gravity, the action includes the Einstein-Maxwell sector Sg with fields collectively 
denoted as g, A in eq. (H)) and the matter sector with the fermion fields 5^, eq. (PT|) . The latter is given as (Euclidean 
signature) : 

Stp = J d^XEy'-gsi^ {D + m) ip, (95) 

where D = T^Vm and the covariant derivative is 'Di\j = Om + \i^abM^°''' — *9^a/- We can always scale away the 
spin connection by redefining the spinor field as in eq. (|A6p . Finite density is described by the electrically charged 
black hole with charge Q, that generates the imaginary time component of the vector potential At^ (cq. [17)) . Radial 
profile of the vector potential Atj^ = /x(l — i) (in dimensionless units) ensures a finite chemical potential at the field 
theory boundary — > /i at r ^ oo, where /i = gpQ (in dimensionless units). Integrating out the fermion fields, the 
gravitational partition function can be written schematically as 

Z= det{D{g,,A,)+m)e-^^^S"^'\ (96) 

where Sg is the Euclidean gravitational action at the saddle points g^,,A^,. The determinant describes fluctuations 
about the saddle point solution , A^. and corresponds to 1 /N correction to the large N limit of a dual gauge theory. 
Because the Euclidean gamma matrices are Hermitian by convention (the signature of the metric fixes the hermiticity) , 
we have = and 

~ with i — 1,2,3, so the covariant derivative is anti-Hermitian. Now it remains to 
show that the determinant of this anti-Hermitian differential operator is real and positive definite [32| . 
Using the anti-commutation relations {r|;,r^*} = 0, we have 



T%Dr% = -D = D\ 



(97) 
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where D = D[g, A). Therefore the determinant 

detD = dct(r|i:>r|) = dctZ?^ = (detZ?)*, (98) 

is real. To show the positive definiteness, we remind the reader that the cigcnmodcs of an anti-Herniitian derivative 
operator come in pairs. If (A,'?/') is an eigenmode of D, 

D^j = XtP, (99) 

then, from eq. (|99|) : 

Dir%i,) ^ {-X)iT%^), (100) 
so (— A,r|;-0) is also an eigenmode of D. Due to anti-hermiticity, from eq. (|97)) : 

i^(r|v) = A*(r|v), (101) 

this eigenvalue is pure imaginary (or zero), — A = A*. The determinant is a product of all the paired eigenvalues, 

det(i:> + m) ^ n,(A, + m){-X, + m) = U^{\X, + mp), (102) 

which is positive definite (or zero). 

In field theory, the eigenmodes of the operator D + fij'^ + m still come in pairs {X,ip) and (— A,7|;^/'). However, 
since is Hermitian, A is no longer pure imaginary, and therefore det(D + ^'^j^+m) is not necessarily positive. The 
sign problem occurs when det is negative for some gauge configurations, or, in other words, it is generically present 
when considering interacting matter at finite density. 



VIII. CONCLUSIONS 



We have studied strongly coupled electron systems in the magnetic field focussing on the Fermi level structure, using 
the AdS / CFT correspondence. These systems are dual to Dirac fermions placed in the background of the electrically 
and magnetically charged AdS-Reissner- Nordstrom black hole. At strong magnetic fields the dual system "lives" near 
the black hole horizon, which substantially modifies the Fermi level structure. As we dial the magnetic field higher, 
the system exhibits the non-Fermi liquid behavior and then crosses back to the conformal regime. In our analysis 
we have concentrated on the the Fermi liquid regime and obtained the dependence of the Fermi momentum kp and 
Fermi velocity vp on the magnetic field. Remarkably, kp exhibits the square root behavior, with vp staying close to 
the speed of light in a wide range of magnetic fields, while it rapidly vanishes at a critical magnetic field which is 
relatively high. Such behavior indicates that the system may have a phase transition. 

The magnetic system can be rescaled to a zero-field configuration which is thermodynamically equivalent to the 
original one. This simple result can actually be seen already at the level of field theory: the additional scale brought 
about by the magnetic field does not show up in thermodynamic quantities meaning, in particular, that the behavior 
in the vicinity of quantum critical points is expected to remain largely uninfluenced by the magnetic field, retaining 
its conformal invariance. In the light of current condensed matter knowledge, this is surprising and might in fact be 
a good opportunity to test the applicability of the probe limit in the real world: if this behavior is not seen, this 
suggests that one has to include the backreaction to metric to arrive at a realistic description. 

In the field theory frame, we have calculated the DC conductivity using kp and vp values extracted from holography. 
The holographic calculation of conductivity that takes into account the fermions corresponds to the corrections 
of subleading order in 1/A^ in the field theory and is very involved [lH|- As we are not interested in the vertex 
renormalization due to gravity (it does not change the magnetic field dependence of the conductivity), we have 
performed our calculations directly in the field theory with AdS gravity-dressed fermion propagators. Instead of 
controlling the occupancy of the Landau levels by changing the chemical potential (as is usual in non-holographic 
setups), we have controlled the filling of the Landau levels by varying the Fermi energy level through the magnetic 
field. At zero temperature, we have reproduced the integer QHE of the Hall conductivity, which is observed in 
graphene at moderate magnetic fields. While the findings on equilibrium physics (Landau quantization, magnetic 
phase transitions and crossovers) are within expectations and indeed corroborate the meaningfulness of the AdS / CFT 
approach as compared to the well-known facts, the detection of the QHE is somewhat surprising as the spatial 
boundary effects are ignored in our setup. We plan to address this question in further work. 

Interestingly, the AdS geometry produces several Fermi surfaces. Theories where the gravity duals have larger 
fermion charge q posses more Fermi surfaces. We find that, in a multi-Fermi surface picture, the Hall conductivity is 
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quantized in a way reminiscent of fractional QHE. By reducing the magnetic field, new Fermi surfaces open up and the 
quantization of Hall conductivity alternates between two different patterns, corresponding to odd and even number 
of Fermi surfaces. It turns out that odd number of the Fermi surfaces results in IQHE plateaus, while even number of 
surfaces gives new plateaus characteristic for the FQHE. In a multi-Fermi surface picture, the quantum Hall plateaus 
show a less regular pattern that agrees with experiments on thin graphite in strong magnetic field (23 | . In our model it 
happens due to the fact that as one of several Fermi surfaces crosses the Landau level, the Hall conductivity jumps to 
a new plateau. This process is not synchronized between different Fermi surfaces. We associate the average distance 
between the Fermi levels with the energy gap usually arising in the FQHE. 

Notably, the AdS-Reissner-Nordstrom black hole background gives a vanishing Fermi velocity at high magnetic 
fields. It happens at the point when the IR conformal dimension of the corresponding field theory is = i, which is 
the borderline between the Fermi and non-Fermi liquids. Vanishing Fermi velocity was also observed at high enough 
fermion charge [2lj. As in [2l|, it is explained by the red shift on the gravity side, because at strong magnetic fields 
the fermion wavefunction is supported near the black hole horizon modifying substantially the Fermi velocity. In our 
model, vanishing Fermi velocity leads to zero occupancy of the Landau levels by stable quasiparticles that results in 
vanishing regular Fermi liquid contribution to the Hall conductivity and the longitudinal conductivity. The dominant 
contribution to both now comes from the non-Fermi liquid and conformal contributions. We associate such change in 
the behavior of conductivities with a metal-" strange metal" phase transition. Experiments on highly oriented pyrolitic 
graphite support the existence of a finite " offset" magnetic field he at T = where the resistivity qualitatively changes 
its behavior [2^. At T ^ 0, it has been associated with the metal-semiconducting phase transition [2^. It is worthwhile 
to study the temperature dependence of the conductivity in order to understand this phase transition better. 

Finally, we suggest as a possibly interesting extension of the current AdS/CFT methodology to compute the gravity 
dual of the finite density matter in Monte-Carlo lattice simulations. This is possible since the sign problem does not 
arise in the holographic setting of a finite density system. Unlike the conventional field theory setup, finite density 
in holography is introduced through an electrically charged black hole, and does not involve matter fields (this is 
also the reason why holography gives universal results: it does not depend on the expectation values of matter fields 
at the leading order). In the gravity geometry, Dirac fermions are coupled minimally to the electric field via the 
covariant derivative. We have shown that the covariant derivative is anti-Hermitian in Euclidean signature, leading 
to the real and positive definite fermion determinant. This makes it possible to simulate finite density systems on the 
lattice in the AdS-gravity geometry using the curved space-time lattice formulation [s^. The simplest holographic 
setup which describes a finite charge density system includes a local U{1) gauge symmetry. Finite density systems 
with global U{1) symmetry can not be simulated numerically in field theory due to the problem with the Gauss law 
in the lattice formulation. Another important advantage of performing Monte-Carlo simulation is that it includes 
the quantum fluctuations for the gauge and gravitational field. So far most calculations have been done in the probe 
limit, with the frozen background for the metric and gauge fields. Analytic calculations which include backreaction are 
usually involved and are done in the next to leading order, e.g. @. Holographic lattice calculations allow to consider 
dynamical gauge and gravity fields with matter, which mimics complicated strong interactions in finite density systems 
and opens a way toward studying novel state of matter and instability mechanisms. 
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Appendix A: Dirac equation in magnetic field 

Here we solve analytically the part of the Dirac equation which depends on magnetic field and space-time coordinates 
of the boundary theory. The free spinor action in the geometry given by eq. ^ and in the presence of magnetic field 
© is given by eq. (PT|) . 

Using the translational invariance, 



tpit, X, y,r)^ / dLjdkc-"^*+'''y k, x, r), 



(Al) 



with k — ky, the Dirac equation (eq. I23p can be written as 



\/~9tt ^ \/9rr \j9ii ^ 



+ T^dk + -UyfyT^^ — iqAy{x)) — m ) ip{u}, fc, x, r) = 0, 



1 



(A2) 



where ga = gxx = 5yy, and At{r) = — ro/r), Ay{x) = hx. From the torsion-free condition, lu^ A = -de", we 
find the spin connection [s^] for the metric [2l 



= — at, ujif. = ^^^dx 



/ 9rr 



/9r 



where i = x,y. Note that 



and 



(A3) 



(A4) 



1 



4'" ''''''' 



1 



/9r 



-r drhi { --^ 

9rr 



rr 4 ■\J~Qr- 

1/4 



where g is the determinant of the metric. Therefore, we can rcscale the spinor field: 

-1/4 



_9_ 

ffrr 



(A5) 



(A6) 



and remove the spin connection completely. The new covariant derivative docs not contain the spin connection so 

^'m = - ^qAm- 

In new field variables, the Dirac equation is given by 



/9v. 



+ i{k-qhx))'^>{uj,k,x,r)=0, 
with = fiq. We separate the x- and r-dependent parts: 



(A7) 



P{r) 



/grr V~9tt 

Q{x) = r^dx+T" i{k-qhx) 



9ii T-.i ■/ , |■^ ^0 NN , 

r l(uj + ^q(l )) - y/Cliim, 

r 



(A8) 



and the Dirac equation is 



(P(r)+g(x))$ = 0. 



(A9) 
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Even though [P{r), Q{x)] 7^ 0, one can find a transformation matrix U such that [UP{r), UQ{x)] = 0, and then look 
for common eigenvectors of UP{r) and UQ{x) as they are commuting Hermitian operators, i. e., the Dirac equation 
reads 



(AlO) 



where I labels the Landau levels. We use / for the Landau index, so as not to confuse it with the Matsubara frequency 
index n. Transformation matrix U should satisfy the conditions 



{u, = 0, {[/, r*} = 0, [u, r*] = 0, [[/, r^] 

which do not fix U completely. It is convenient to use the following basis Q: 



0, 



(All) 



— (7 

- 

a2 
a' 







F5 = 



ia^ 

i(7^ 

ia^ 
-ia^ 



F^ = 







Note, that the following relation holds 
as expected, with ^ t, 1 ^ x, 2 ^ y, 3 



p5 j-^0-plp2-p3 

r. In the representation of eq. (jA12[) . we can choose 




U 



-la 




(A12) 
(A13) 
(AM) 



We split the 4-component spinors into two 2-component spinors (we do not write zero entries) F = (Fx^F^)^ where 
the index a = 1, 2 is the Dirac index of the boundary theory, using projectors 



n„ = -(i-(-i)"F'-F^F^), a = 1,2, ni + n2 = i. 



(A15) 

which commute with the Dirac operator of eq. ([37| . and Fa = IIq.^, a = 1,2, decouple from each other. Gamma 
matrices in eq. (jA12p were chosen in such a way that this decoupling is possible. 
Writing the Dirac equation (jAlOp for F = (Fi, F^)^ , we have 



'guicr m 



idx + A; (fc — qhx) 
{k - qhx) -idx + Xi 



F2 



cr^(w + fj.q{l - ro/r)) - A„ 



= 0, 



Fi 
F2 



= 



(A16) 



where in X (E)Y, X acts inside Fi or F2 and Y acts between Fi and i^2- In eq. (jA16|) . the 1 in the first equation 
shows that there is no mixing of Fi and F2 by the operator UP{r) and the 1 in the second equation shows that there 
is no mixing inside Fi or F2 by the operator UQ{x). Therefore, the solution can be represented as 



F2 



f fi'\r)gl'\x)\ 
fl^\r)9['\x) 
fl'\r)9f\x) 

\fl^\r)gf\x) J 



(A17) 



We do not write explicitly the dependence on w and k. It is convenient to make a unitary transformation: 



= M 



.(2) 



M = 



1 -i 
-i 1 



(A18) 



Dirac equations for each component are written as: 



/9i. 



9rr 
V9i: 



l9v. 



\J-9tt 
V~9tt 



iio + fi,{l-ro/r)) + Xi],f\r) = 0, 



ioj + t,,il-ro/r)) + Xijf^'\r) = 0, 
(a,+5;)^(2)_A;C(i) = 0. 



(A19) 
(A20) 
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In equations IA20I for the x-dependent part, we have rescaled x - 
second order ordinary differential equations 



\qh\ {x — ^) with k = ky and A; = ^/\qh\ A;. The 



(A21) 



witli p — 1,2, arc solved by substitution C^''^ = e~^ /2^(p), This is exactly the Schrodinger equation for a harmonic 
oscillator, so the eigenfunctions are Hermite polynomials and we obtain the following solutions, indexed by an integer 
I e Z that is related to the eigenvalue A; by A; = y/2\qh\l: 



(A22) 

The normalization constant Ni is proportional to l/VWll. Substituting the solutions from eq. (|A22p into the first 
order eigenvalue equation with x-dependence gives the following solutions: 



F, 



and 



f\r)(:l'\x) 



Xl = V^\qh\i, 



(A23) 



F, = 



f?\r)cl'\x) 
-fl'\r)cl'\x) 



A; = -^J2\qh\^. 



(A24) 



Solving the first order x-depcndcnt equation, we get the same eigenvalue, but slightly different eigenfunctions for 
different signs of qh. In particular, e. g., for F, the pairs /^^•'(C^^^ C*'^^)'^ ^^nd /'■^•'(C''^\ — C'"^'')"^ correspond to g/i > 
and qh < 0, respectively Different sign of qh stands for the positive and negative Landau level index I. 

Finally, the general solution to the Dirac equation is given by a linear combination of eqs. (|A23p and (|A24|) : 



F,,i=J2iaiFi+biFi). 
I 

Using the eigenvalues determined by eqs. (|A23IIA24|) in the equation for the radial part (|A19|) . we get 



(A25) 



z<y'^dr — m 



\J-9tt 



a^{uj + ^q(l - ro/r)) 



'yj2\qh\l \ ® 1 



F2 



= 0, 



with I ~ 0, l,dots; and the same for F replacing y^2\qh\l —y^2\qh\l. It coincides with & 
equation at zero magnetic field) with the momentum replaced by the Landau level eigenvalue [2: 



(A26) 
(A14) in (Dirac 



k ±^/2\qh\l. (A27) 
Equation (|A27p also gives a prescription on how to treat the limit of zero magnetic field — > 0. The limit is to be 



taken keeping, e. g., 2\qh\{l + 1) = kp fixed as /i — >■ 0. In a compact form, the Dirac equation in a magnetic field (|A7j 
is given by 



1 



1 



F* i{u: + qAt) - m 



-.u-'yAmnF{r) = o, 



(A28) 



with F = {F i,F2Y , I 
given by eq. (jA14P : 



0, 1, . . ., for F replace ^2\qh\l 



u- 



-^2\qh\l, and J7 ^ is the matrix inverse to the matrix 

(A29) 



which we use in the main text. 
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Appendix B: Spectral function 

111 what follows we use the dimensionless variables (|15m7|) . Following the analysis of Q, the flow of the Green's 
function is determined by 



Gi?(a;,0 = lime 









(Bl) 



where C+ ('') = ^oJ and (r) = from the solutions (|A23IIA24| . In obtaining this relation, we absorbed the 

coefficients appearingin eq. (|A25|) into the definitions of the radial functions. The functions ^j^-* satisfy the following 
differential equation Q: 

^drd^ = -2mV^d" + {u{r) ± Xind^ f + (^M T A,), (B2) 




Qrr 

with u{r) given by 



uir) = J^iu + qAtir)). (B3) 
V -9tt 

Writing explicitly in the metric eq. (fT8|) . we have 

T^^/ldrSH^ = -2mr^l^ + {u{r) ± A/)(e^V + Hr) T A,), (B4) 

where u(r) is given by 

«(r) = -^(c. + ^,(I-i)) (B5) 
with / = (''-1) ('■^+2r+3) Q^j^ — Q^ Near the horizon (r = 1) the flow equation reduces to 



r29.e^) = i(d'^ + lf, (B6) 



which due to the double zero in / has a regular solution only if ^±(r = 1) = ±i. Writing the radial equation in terms 

of ^ and choosing the infalling boundary conditions fixes £!± {r = 1) = i. 

The key quantity that we extract from the Green's function is the fermionic spectral function 

A{uj, ky) = Tr (ImGfl(cj, ky)) , (B7) 

which we analyze in the main text of the paper. 
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